Magnetic soft modes in the locally distorted triangular antiferromagnet a-CaCr 2 04 
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In this paper we explore the phase diagram and excitations of a distorted triangular lattice 
antiferromagnet. The unique two-dimensional distortion considered here is very different from the 
'isosceles'-type distortion that has been extensively investigated. We show that it is able to stabilize 
a 120° spin structure for a large range of exchange interaction values, while new structures are found 
for extreme distortions. A physical realization of this model is a-CaCr2 04 which has 120° structure 
but lies very close to the phase boundary. This is verified by inelastic neutron scattering which 
reveals unusual roton-like minima at reciprocal space points different from those corresponding to 
the magnetic order. 



PACS numbers: 61.05.F-, 75.10.Hk, 75.30.Ds, 75.30.Et 

In frustrated antiferromagnets (AF) the magnetic in- 
teractions compete and it is impossible to satisfy them 
all simultaneously. This competition causes degenerate 
ground states to occur, and the system fluctuates be- 
tween them suppressing long-range magnetic order to 
temperatures well below the Curie- Weiss temperature [1] . 
Most frustrated magnets do eventually order and in the 
case of isotropic or Heisenberg interactions a compromise 
is often reached where the magnetic interactions are all 
partially satisfied by a non-collinear arrangement of the 
magnetic moments. 

The triangular lattice Heisenberg antiferromagnet 
(TLHAF) is one of the simplest examples of an extended 
frustrated magnetic system, however open questions re- 
main concerning its physical properties. In the case of 
spin-1/2, the ground state is expected to show spin-liquid 
behavior and possibly form a resonating valance bond 
(RVB) state [2]. For spin-1, both a RVB state and a spin 
nematic phase have been proposed [2j [3] . Systems with 
larger spin values order in a helical structure character- 
ized by a magnetic ordering wavevector of 1/3 where the 
angle between nearest neighbor spins is 120°. 

Linear spin-wave theory (LSWT) for the TLHAF pre- 
dicts one sinusoidal, doubly degenerate mode. How- 
ever the excitation spectrum for spin-1/2 is expected to 
strongly deviate from LSWT; downward renormalization 
(of up to 33%) and damping of excitations are expected 
over large regions of the Brillouin-zone [3H6]. A reduced 
effect is expected for higher spin values. Introduction of 
antiferromagnetic next-nearest-neighbor interactions also 
modifies the excitations and drives the magnetic struc- 
ture from helical to collinear when the ratio of next- 



nearest to nearest neighbor interactions is greater than 
the critical value of 0.125 [7 j. In the collinear phase, 
quantum fluctuations lift the classical degeneracy and se- 
lect collinear stripe order via the mechanism of order by 
disorder [8 , they also significantly alter the spin- wave 
excitation spectra around certain regions of the Brillouin 
zone [7]. 

Most physical realizations of triangular lattices are dis- 
torted. The distortion usually constitutes a shortening 
of the nearest neighbor distances along one of the three 
equivalent directions, resulting in one of the exchange 
constants being larger than the other two. This distor- 
tion partially lifts the frustration and promotes the for- 
mation of chains, with frustrated interchain coupling as 
in Cs 2 CuCl 4 [9 and CuCr0 2 [10 . The resulting order- 
ing wavevector is intermediate between the undistorted 
value of 1/3 and the AF chain value of 1/2. 

In this paper, we investigate the excitation spectrum of 
a-CaC^CU which is another distorted triangular lattice 
antiferromagnet but where the distortion is very different 
from those investigated previously. Furthermore, in ad- 
dition to interactions between nearest neighbor magnetic 
ions there are also substantial next-nearest-neighbor in- 
teractions. This system allows examination of the effects 
of competing nearest- neighbor and next-nearest-neighbor 
interactions on a distorted triangular lattice where the 
distortion places the system close to a magnetic phase 
boundary. 

a-CaCr 2 04 belongs to the orthorhombic Pmmn space 
group and has a distorted delafossite structure [11 . The 
magnetic Cr 3+ ions lie in an octahedral environment 
and are characterized by a half-filled t 2g shell giving 
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FIG. 1. (color online) (a-b) Ordered magnetic structure of 
single triangular layer of Cr 3+ ions, arrows represent the mag- 
netic moments and spins are rotated to the be plane for better 
visibility, (a) and (b) give that structure at the position of the 
square (km = (0,1/3,0)) and circle (k m = (0,0,1/2)) shown 
on the phase diagram (d) respectively. The solid lilac and 
yellow lines represent J zz i and J ZZ 2, the dashed lines are J c hi 
and J c h2 respectively. Thin gray lines depict next-nearest- 
neighbor interactions, (c-e) Magnetic ordering wavevector as 
a function of Ai, A2 and J nnn for A3 = 0, the colors indicate 
the size of the components of k m , in the shaded region the 
magnetic ground state has a multi-A; structure. 



rise to spin-3/2 and quenched orbital angular momentum 
which ensure that the magnetic interactions are isotropic. 
They lie at two inequivalent crystallographic positions 
and form distorted triangular layers stacked along the 
a axis. There are two layers and eight magnetic ions 
per unit cell. Together with a-SrC^C^ [T2] , this com- 
pound shows a unique distortion, see Fig.ga). Although 
the Cr 3+ ions are locally displaced from the ideal trian- 
gular positions, the ratio of the b and c lattice param- 
eters remain close to the triangular 2/y/3 value. The 
four inequivalent nearest neighbor Cr 3+ -Cr 3+ distances 
range from 2.889 A to 2.939 A and form two zig-zag and 
two straight chains. The short Cr 3+ -Cr 3+ bonds and 
the almost 90° Cr 3+ -0 2_ -Cr 3+ angles mean that the 
dominant exchange mechanism is direct exchange as in 



many other Cr 3+ delafossite compounds [TDJ [T3HT9] . The 
Curie- Weiss temperature of -564(4) K, yields an average 
nearest neighbor exchange constant of J nn = 6.48 meV 
[20] [21]. However the sensitive distance dependence of di- 
rect exchange would suggest large differences between the 
inequivalent nearest-neighbor interaction strengths. The 
interplane distances are much greater (~5.524 A) and 
are therefore expected to be mediated by much weaker 
superexchange interactions. 

a-CaC^C^ orders magnetically below 7n = 42.6 K 
[20] |2T] . The ordered state is the same as for the 
undistorted TLHAF with angles of 120° between near- 
est neighbors and spins lying in the ac plane, (Fig. [TJa)). 
In ref. 21, a classical model was developed which showed 
that this 120° structure is stable over a large range of 
the exchange constant values and provided constraints on 
their relative sizes. In this paper, we show that despite 
the simplicity of the magnetic structure, the magnetic 
excitations of a-CaCr204 are highly complex. Further- 
more, a significant softening of the spin- waves was found 
around wavevectors such as (0,1/2,0) and (0,0,1) placing 
it close to an instability away from the simple 120° struc- 
ture towards the formation of a complex multi-/c struc- 
ture. 

Both polycrystalline and single-crystal samples of a- 
CaC^CU were measured [22]. The polycrystalline sample 
has a mass of 3.5 g; the single-crystal is 340 mg and 
consists of three crystallographic twins rotated by 60° 
with respect to each other about the shared a-axis. 
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FIG. 2. (color online) (a) Energy and wavevector dependence 
of ARCS spectra for E\ = 40 meV. Right side shows an energy 
cut for 1.5<Q<3.3A from the Ei = 60 meV data measured 
at 4.2 K (black line) and 55 K (orange line). Arrows indicate 
the van Hove singularities at 4.2 K. (b) Powder average of 
calculated spectra for a spatially isotropic triangular lattice 
with interplane coupling, dashed line indicates the border of 
the measured region. The right side shows the equivalent 
energy cut. Intensities are in arbitrary units. 

Powder inelastic neutron scattering was measured on 
the ARCS [23] chopper spectrometer at SNS, ORNL, 
USA. Data was collected with incident neutron energies 
of 40 meV and 60 meV (resolution of 1.64 meV and 2.35 
meV at the elastic line respectively) at 4.2 K and 55 K. 
To map out the excitations along all reciprocal space di- 
rections, a single crystal experiment was performed on 
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the Merlin [24] chopper spectrometer at the ISIS facility, 
RAL, UK. The sample was measured at 5 K using an 
incident neutron energy of 50 meV. The crystal was ro- 
tated over 120° in 1° steps about the vertical a axis and 
data was collected for 1 hour per angle. The data was 
transformed into units of energy and wavevector trans- 
fer, combined and integrated over h using the Horace 
software [25 . Single crystal experiments were also per- 
formed on the Puma thermal triple axis spectrometer at 
FRM2, Germany. The final neutron wavelength was fixed 
to 2.662 A -1 , giving an energy resolution of 0.75 meV. 
The magnetic excitations were measured at 3.5 K, while 
several scans were repeated at 220 K to separate phonon 
and magnon scattering. To determine the strength of 
the weak interplane coupling and check for an anisotropy 
gap, scans were performed on the cold neutron high res- 
olution V2 triple axis spectrometer at HZB, Germany. A 
fixed final neutron wavelength of 1.20 A -1 gave an energy 
resolution of 0.093 meV. 

The powder spectra shows strong inelastic inten- 
sity centered around Q~1.54 A -1 corresponding to the 
(1,4/3,0) magnetic Bragg-peak (1.55 A -1 ), see Fig.^a), 
where the characteristic V-shaped dispersion is visible. 
Van Hove singularities are observed at 5, 11, 17 and 33 
meV indicating energies where the dispersion is flat, how- 
ever no sharp upper edge of the dispersion is observed 
as would be naively expected in the classical case. The 
data can be compared to a spin wave simulation, as- 
suming spatially isotropic nearest-neighbor interactions 
J nn =6.48 meV (from the high temperature susceptibil- 
ity) and interplane coupling J[ nt =0.005 meV. The sim- 
ulated spectra in Fig. [2^b) shows van Hove singularities 
at 19 meV and at the top of the magnon dispersion at 
32 meV. These singularities are missing in the measured 
spectra and the magnon intensity decreases smoothly to 
background with increasing energy transfer. It is clear 
that despite its 120° magnetic structure the excitations 
of a-CaC^C^ cannot be approximated by an isotropic 
triangular antiferromagnet and the spatial anisotropy of 
the interactions must be considered. 

The full complexity of the spectra is revealed by the 
single crystal data. Figure [3] (a) shows the (0,4/3,/) 
dispersion from Merlin. The scattering intensity is 
strong around the magnetic zone centers (0,4/3, ±2) and 
(0,4/3,0) as expected. Many spin- wave branches dis- 
perse from these points, the steepest of which reach a 
zone boundary energy of 35 meV. Details of the lower 
energy branches are revealed by Puma. Along (0,k,l) 
(Fig. [3jb)), the strong scattering around the zone center 
is again visible at (0,2/3,1), the V-shaped cone is typical 
of the magnon dispersion of the spatially isotropic trian- 
gular lattice. However it is broader than the resolution 
probably due to multiple spin- wave branches arising from 
having eight magnetic ions in the unit cell and three crys- 
tallographic twins. The gap at (0,2/3,1) is due to the in- 
terplane coupling. On both sides of the cone, low energy 
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FIG. 3. (color online) (a) Spectra measured on Merlin, (b- 
c) spectra measured on Puma, black dots denote measure- 
ment points, (d-f) Calculated spin- wave spectra using the 
best fit parameters (Tab. [I]). The instrumental resolution is 
convolved with the calculated spectra, using the resolution 
ellipsoid obtained from Rescal [26] (Popovici's method [27]). 
To reproduce all the observed modes, the dispersion of the 
three crystallographic twins were combined, weighted by their 
experimentally determined volume fractions. Intensities are 
in arbitrary units, , "r.l.u." stands for reciprocal lattice units. 



modes (E < 10) are visible with local minima at (0,0,1) 
and (0,4/3,1). The scattering in the region 1.3 < k < 2 
and 15 < E < 25 is from phonons as confirmed by com- 
paring to the spectra at high temperatures (T=220 K). 
The low energy modes are also visible along (0,0,/), (Fig. 
[3] (c)). They are especially strong around / « 1 and have 
a distinctive shape. The broad modes above 15 meV are 
again phonons. 

In order to probe the interplane coupling, energy scans 
at the zone center were repeated at different values of out- 
of-plane wavevector. The scan at the (1,4/3,0) magnetic 
Bragg peak shows no gap revealing that anisotropy is 
negligible. However the magnetic excitations at (2,4/3,0) 
have a gap of 3.5 meV revealing a small but significant 
antiferromagnetic coupling between adjacent planes. 

A spin-wave model was developed to explain the ob- 
served features in the excitation spectra of a-CaC^C^. 
Based on the crystal structure and magnetization data, 
the magnetic properties should be well accounted for by 
a Heisenberg model. Four inequivalent nearest neighbor 
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and four inequivalent next-nearest-neighbor interactions 
in the triangular plane are allowed by symmetry. The 
interplane exchange interactions are assumed to be all 
equal. We consider the spin Hamiltonian to be com- 
posed of interplane, nearest-neighbor and next-nearest- 
neighbor terms. Since the system is strongly frustrated, 
the magnetic ordering wavevector and spin directions are 
sensitive to the sizes of the interactions, for details see 
Ref. 21. The observed value of the ordering wavevec- 
tor fc m =(0, 1/3,0) gives the following restriction on the 
parameters: 

A3 = J zz i + J zz2 — J c hl — Jch2 = 0, (1) 

where the zz and ch indices denote the zig-zag and chain 
nearest-neighbor in-plane interactions respectively. The 
next-nearest-neighbor interactions ( J n nni-^nnn4) do not 
influence the ordering wavevector. The above constraint 
is applied in the rest of this paper. We define a set of 
parameters linearly independent of A3: 

Jnn ~ ^(^chl + J c h2 + ^zzl + ^zz2)? 

Ai = — — (J zz i — J zz2 ), A 2 = — - — (J c hl — Jch2)- (2) 
ZJ nn ^nn 

The magnetic structure as a function of the new parame- 
ters is determined by minimizing the energy, assuming a 
single pair of modulation wave vectors k m and — k m [28] . 

Figure [ljc-e) shows the phase diagram as a function 
of Ai, A2 and average next-nearest-neighbor interaction 
J nnn . The 120° structure characterized by fc m =(0, 1/3,0) 
is stable over a large range of exchange constants. In the 
shaded region the magnetic ground state has a multi-A: 
structure, which is beyond our current model. A Monte 
Carlo method applied to an extended unit cell would be 
necessary to determine the structure here. Beyond the 
shaded region the magnetic structure can again be de- 
scribed by a single wavevector but with a new value e.g. 
fc m — (0, 0,1/2) represented by the circle in Fig. [TJd) and 
illustrated in Fig. [TJb) 



TABLE I. Values of the fitted exchange parameters. 



Name 


Jzzi 


Jzz2 


«/<chl 




«/int 


Value (meV) 


5.8(6) 


11.8(14) 


9.1(10) 


8.6(15) 


0.027(1) 


Cr-Cr dist. (A) 


2.939 


2.889 


2.911 


2.907 


5.529 


Name 


t/nnnl 


t/nnn2 


«/nnn3 


<^nnn4 




Value (meV) 


0.11 


0.31 


0.96 


0.90 





Linear spin- wave theory was used to calculate the exci- 
tation spectra and neutron scattering cross section of the 
proposed Hamiltonian [6j[29j[30]. LSWT fails if the mag- 
netic ground state structure is incorrect for a given set 
of exchange constants, and since a single pair of ordering 
wavevectors is assumed in the spin-wave calculation, the 



borders of the multi-A: phase can be determined by the 
values of the exchange interactions where the spin-wave 
energies become imaginary, (shaded regions in Fig. [TJc- 
e). Using the parameter constraint A3 = 0, we fitted 

the parameters J nn , Ai, A 2 , J in t, Jnnnl, ^nnn2, */nnn3, 

Jnnn4- The calculated modes with strongest intensity 
were fitted to the dispersion curves extracted from the 
data. The fit favored a set of exchange interactions close 
to the boundary between the k m = (0,1/3,0) and multi- 
k phase. In fact the fit was unstable, driving the ex- 
change parameters out of the 120° region, and it had to 
be constrained to stay within this phase to ensure con- 
sistency with the diffraction results. The best fits were 
obtained with the following values: J nn = 8.8(8) meV, 
Ai = ±0.34(4), A 2 = 0.03(6) and the average next- 
nearest-neighbor strength J nnn = 0.69 meV. Since the 
four inequivalent next-nearest-neighbor interactions are 
strongly correlated, no unique set of values were found. 
A good solution is shown in Tab. [I] and the corresponding 
fitted dispersion relations are shown alongside the data 
in Fig. [3] 

The values of Ai and A 2 correlate well with the Cr 3+ - 
Cr 3+ distances determined by neutron powder diffrac- 
tion. The larger Ai value arises from the difference be- 
tween the two zig-zag interaction strengths, where the 
difference between their distances is 0.050 A. In con- 
trast, the difference between the two chain distances is 
only 0.004 A, thus explaining the very small difference 
between their interactions (A 2 ). The sign of Ai is unde- 
termined, because it does not affect the spin-wave spec- 
tra, however to make the shorter coupling stronger Ai 
has been set negative. These values for the first neighbor 
exchange constants agree well with the Cr 3+ -Cr 3+ in- 
teraction strength as a function of distance as found for 
other chromium oxides delafossites with direct exchange 
interactions [T0l[T3Ul9]. 

The results of the LSWT fit place a-CaCr 2 4 at the 
edge of the 120° phase close to the multi-/c phase, as 
shown by the circle in Fig. [ljd) . The proximity of the 
phase boundary is revealed by the presence of low energy 
modes with roton-like minima at wavevectors that do not 
correspond to the k m = (0,1/3,0) phase but rather act 
as soft modes of the nearby structure. LSWT is able to 
reproduce the spectra but there are small discrepancies 
(Fig. [3]). It has been shown theoretically that for the 
isotropic triangular AF the non-collinearity of the spin 
moments gives rise to substantial downwards renormal- 
ization and broadening compared to LSWT and while 
these effects are most pronounced for spin- 1/2 [4]-[6], re- 
cent calculations suggest that they may be ~10 % for 
spin-3/2 [31]. Quantum effects are also enhanced for low 
energy modes close to a phase boundary. Soft modes are 
predicted for the collinear phase of the isotropic triangu- 
lar lattice with next-nearest-neighbor interactions [7, 8 
corresponding to several possible phases which are classi- 
cally degenerate. Via the mechanism of order by disorder 
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quantum fluctuations (QF) select a striped collinear state 
by keeping the relevant soft mode gapless while increasing 
the gap of the other modes. It was also shown, that QF 
can shift the magnetic phase boundary of the isotropic 
triangular AF with next-nearest-neighbor interactions in 
external magnetic field [32] . 

Since QF have a significant effect on the isotropic 
triangular AF close to phase boundaries, it is possible 
that they also modify the spectra of a-CaCr204. These 
modifications might be responsible for stabilizing the 
120° phase via the mechanism of order by disorder even 
though the LSWT fit to the dispersion suggests that 
a-CaCr204 might lie outside this phase. It should be 
noted that the isostructural compound a-SrCr204 also 
has a 120° magnetic structure, here the crystal structure 
suggests that the difference between chain interactions is 
greater than in a-CaC^CU while the difference between 
the zig-zags in smaller [12]. The excitation spectra of 
this compound has not been measured. Further theoret- 
ical study taking quantum fluctuations into account is 
necessary to determine the phase boundary and excita- 
tion spectra. 

In conclusion, despite its highly symmetric 120° mag- 
netic structure, the magnetic excitation spectrum of a- 
CaCr204 is highly complex. The unique pattern of near- 
est neighbor exchange interactions as well as the substan- 
tial next-nearest-neighbor interactions place it close to 
the phase boundary of the 120° structure as is clearly re- 
vealed by the presence of low energy modes acting as soft- 
modes of the neighboring structure. Indeed the LSWT 
fit tends to drive the system towards the nearby mufti- 
phase in contradiction to the observed ordering wavevec- 
tor of k Yn = (0,1/3,0) and quantum fluctuations may be 
necessary to stabilize the system at a well-defined point 
inside the 120° phase. 
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